Probing multipartite entanglement in a coupled Jaynes-Cummings system 
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We show how to probe multipartite entanglement in N coupled Jaynes-Cummings cells where the degrees 
of freedom are the electronic energies of each of the N atoms in separate single-mode cavities plus the N 
single-mode fields themselves. Specifically we propose probing the combined system as though it is a dielectric 
medium. The spectral properties and transition rates directly reveal multipartite entanglement signatures. It is 
found that the Hilbert space of the N cell system can be confined to the totally symmetric subspace of two states 
only that are maximally-entangled W states with 27V degrees of freedom. 
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I. INTRODUCTION 

A single two-level atom coupled to a single cavity, or res- 
onator, mode has been studied intensively since the introduc- 
tion of this 'Jaynes-Cummings (JC) Model' independently by 
Jaynes and Cummings jl] and by Paul J2l in 1963. Quantum- 
field effects such as periodic spontaneous collapse and re- 
vival are now studied and observed in many systems and 
as a multitude of manifestations J3l ■ Coupled JC systems have 
been proposed as a basis for quantum networks Q5|] and could 
behave as a novel condensed-matter system if enough JC sys- 
tems can be coupled together J6| . 

In a preliminary study we developed theoretical tools for 
calculating the spectrum, stationary states and dielectric sus- 
ceptibility Q. Here we use and extend those tools, espe- 
cially to include the nontrivial open-system effects of spon- 
taneous emission and cavity losses. In particular, we consider 
the problem of how to probe and characterize such systems 
experimentally. Two quite different approaches are evident. 
One approach is to probe each component in microscopic 
paradigm, namely drive and detect the various atoms and cav- 
ity modes. Another approach, which we favour, is to treat the 
coupled JC system as a 'black-box' model and probe it as a 
single unit following a macroscopic paradigm. 

Our concept is to regard the coupled JC system as a dielec- 
tric medium whose susceptibility carries a signature of the pe- 
culiarities of the coupled JC system. Diagonalizing the Hamil- 
tonian of the coupled JC system, we find entangled states and 
their energies. Using the Fermi's golden rule, we calculate 
the transition rates between different manifolds of the energy 
states of the system. In this paper, we restrict the calculation 
to transitions from the single-excitation states to the ground 
state of the system. In particular we show that the dielectric 
susceptibility of the coupled- JC medium reveals, by probe- 
field spectroscopy, quadripartite entanglement of the system 
comprising mutually coupled atoms and cavity modes. Based 
on our theoretical framework for JC systems mutually coupled 
by overlapping extra-cavity longitudinal fields, we can calcu- 
late stationary states for the coupled system, energy spectrum, 
and dielectric susceptibility. 

The paper is organized as follows. In Sec. [XT] we give a 



qualitative discussion and a detailed calculation of radiative 
properties of a single JC cell. Section|III]is devoted to the dis- 
cussion of the entangled and radiative properties of two cou- 
pled JC cells. We derive single-excitation states of the system 
and show that they are the W state class corresponding to a 
superposition state of single excitations amongst each of the 
two atoms and two modes. We then find under which condi- 
tions the states reduce to the maximally entangled four-qubit 
W states, and how to quantify the degree of entanglement of 
the states using probe-field spectroscopy methods. A general- 
ization of the calculation to an arbitrary number of mutually 
coupled JC cells is presented in Sec. [TV] Finally, in Sec.|V] we 
summarize our results. 



II. ENTANGLEMENT AND RADIATIVE PROPERTIES OF 
THE SINGLE JC CELL 

The JC cell (as we refer to a single JC system) is a com- 
posite atom cavity-mode system that radiates via two dis- 
tinct channels. The closed coupled JC system is character- 
ized by just two parameters, namely the atom-cavity coupling 
rate g and the inter-cavity coherent hopping rate k. One radia- 
tive channel arises due to coupling of the atom to the cavity 
side modes (free-space modes) causing spontaneous decay of 
the atomic excitation with an inhibited spontaneous emission 
rate 7 a . The inhibition refers to the fact that some free radia- 
tive modes are suppressed by cavity confinement of the atom. 
The other radiative decay mode is due to losses from one or 
both mirrors of the optical cavity at rate j c . We assume inde- 
pendence between these two radiative channels. 

A single JC cell comprises a single cavity field mode of 
frequency uj c and a two-level atom with energy states \gi) 
and \e.i) and corresponding energy difference 
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The atom is coupled to the cavity mode with coupling con- 
stant g, which we choose to be real with no loss of general- 
ity. We assume that the atomic transition frequency is detuned 
from the cavity frequency by A = w a — cj c . 

In the rotating-wave approximation, and using a system of 
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units in which H = 1, the Hamiltonian of the single JC sys- 
tem is 

H 3C = uj c yd)a + + icj a (T z + ,g (a^o-_ + aa + ) , (2) 

with and a creation and annihilation operators of the cavity 
mode, respectively, and a± and a z , with spec(<7*) = ±1, the 
spin operators for the atom. 



A. Energy spectrum 

The JC energy spectrum comprises a ground state of en- 
ergy ojq = — A/2 and a ladder of doublets of energies 



uj± n = nuj c ± ng 2 + -A 2 , n=l,2,... (3) 

The corresponding stationary energy states are the ground 
state \0g) and the excited state doublets 

|+, n) = cos 0„ |e, n — 1) + i sin0„ \g, n) , 

|— , n) = — sin 0„ |e, n — 1} + i cos9 n \g, n) , (4) 

with 

20 n = tan"* (^) . (5) 

Note that the doublets are in a pure, bipartite entangled state 
of the atom and the field. This entanglement is maximal (for 
all the doublets) only when A = 0. 



(a) (b) 




-10 -5 5 10 -10 -5 5 10 



A [arb. units] A [arb. units] 

FIG. 1: (Color online) Variation of the energy difference 5uj±i — 
u±i — Ld c of the n = 1 eigenstates with the detuning A (in arbitrary 
units); 8oj-\ (black solid line) and 8lj+i (red dashed line) for a single 
JC cell with: (a) g = and (b) g = 1. 

Figure[TJshows a variation of the energies oj± n of the n = 1 
doublet |±, 1) with the detuning A around the unperturbed 
energy oj c . Notice the level crossing in the absence of the cou- 
pling g at A = 0, and the appearance of the familiar avoided 
crossing effect when g ^ 0. 

In the following we explore the radiative properties of this 
single JC cell, in particular looking for signatures of this atom- 
field entanglement in the transition rates between states with 
different total number of excitations. Specifically, we search 
for conditions that reveal whether the entanglement of the 
eigenstates (O is maximal or not. 



B. Transition rates 

As we have seen, each doublet has some fixed number of 
quanta shared between the atom's electronic state and the cav- 
ity mode field state. We employ the integer v to designate the 
number of quanta in the system. The ground state corresponds 
to v = 0, which means that the electronic state and field state 
are both in the lowest level. The first doublet has v = 1 quan- 
tum, the second doublet has v — 2 quanta, and so on. 

Consider first the case n = 1 with transitions from the 
single-excitation states |±, 1), with splitting due to the atom- 
field interaction, to the ground state, which is a product state 
not affected by the atom-field interaction. The transition rates 
from the excited states to the ground state are given by Fermi's 
golden rule Jl 

r ±il = 7a |(1, ±| a+ \g, 0) | 2 + 7 c | (1, ±| \g, 0) | 2 • (6) 

These transition rates are a sum of transitions caused by spon- 
taneous emission from the atom, with rate j a , and by damp- 
ing of the cavity mode, with rate j c . Consequently the co- 
efficient 7 a quantifies the amount that the atomic spontaneous 
emission contributes to the transition probability. Similarly, j c 
quantifies how much the cavity losses contribute to the transi- 
tion probability. 

From Eq. (|4]i we readily calculate the transition dipole mo- 
ments between the states |±, 1) and \g, 0) whence we obtain 
the transition probabilities 

= 7c + (7a ~ 7c) COS 2 01, 

r-,i= 7a -(7a-7c)cos 2 0i. (7) 

A number of interesting properties are immediately evident 
from the expressions for these probabilities. For example, in 
the absence of either atomic spontaneous emission ( 7a = 0) 
or cavity dissipation ( 7c = 0), that is when only a single dis- 
sipation channel is present in the system, the transition prob- 
abilities depend on the nature of the states from which they 
originate with this nature given by the parameter 0i appear in 
Eqs. © and Q. 

This dependence on 0i causes the two transition proba- 
bilities to be mutually correlated. For example, an increase 
of r + ,i implies a decrease of r_ i and vice versa. The prob- 
abilities become independent of each other and their magni- 
tudes equalize only when the states become maximally entan- 
gled. Thus, we may infer from the transition probabilities to 
what degree the states are entangled. However, this conclu- 
sion is based on a simplified model of the system involving 
only a single decay channel, which is incompatible with what 
one encounters in practice. A practical JC system radiates 
through both dissipation channels. 

We distinguish here two parameter regimes with a qualita- 
tively different behavior for the transition probabilities. These 
two regimes are distinguished according to the magnitude of 
7 C relative to 7a . When the dissipation rates are equal, 7 a = 
7 C = 7, Eq. (0 yields that the probabilities are independent 
of each other and have equal magnitudes, T + ,i = r_ i = 27, 
regardless of the state of the system (i.e. 0i). In other words, 
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the transition probabilities tell us nothing about the nature of 
the states involved. In physical terms, this is a consequence of 
the fact that, for j a — j c , one cannot distinguish which dissi- 
pation channel is used by the quantum leaving the JC system 
regardless of the initial state. Therefore, determining the en- 
tanglement of the states |±, 1) is not possible from transition 
rates if 7 a = -f c so of course the set-up must then avoid this 
condition for entanglement to be measurable. 

For 7 a 7^ 7c, the transition rates depend explicitly on the 
amplitudes of the states involved thereby enabling the de- 
gree of entanglement to be discerned from the transition rates, 
for example by measuring the difference or the ratio between 
r + ,i and r_ t i. From Eq. Q we see that, for non-maximally 
entangled states (cos 2 9\ ^ 1/2), equality between r+.i and 
r_.i cannot be achieved. However, when cos 2 0\ = 1/2, 
which corresponds to the case of maximally entangled states, 
r + .i and r_ i have the same magnitude. Equality between 
the transition probabilities with 7 a ^ j c can only occur for 
maximally entangled states. 

We now consider the case n > 2. In this case, transitions 
occur between two neighboring doublets of entangled states 
(n = 2 and n = 1). Notice that the transitions occur between 
states of different degree of entanglement, and there are two 
possible transition channels from each state of the upper dou- 
blet to states of the doublet below. Transitions from the i th to 
the state of the neighboring n and n — 1 doublets occur 
with rate 



7a|(Ml<7 + \j,n- 1)1 +y c \(n,i\a*\j,n - 1)1 



(8) 



Using Eqs. and (0, we readily find that the transitions oc- 
cur with probabilities 

r+,±n = 2 t(7a _ 7c) T (7a + 7c) cos26> n _i] cos 2 6 n 

+ 2 n ^ c (1 COS ^ n C0S ^n-l) 
1 



± -7c \Jn[n - 1) sin 20 n sin 20 n _i , 
r_,±n = 2 K^a ~ 7c) T (7a + 7c) cos26>„_i] sin 2 6 n 
+ 2 n 7c (1 T cos 26 n cos 20 n _i) 



thereby yielding 

r +i „ = ^ l+J" = n 7c + (7a - 7c) COS 2 9 n , 
J = ± 

r_, n = ^ r_j„ = n^ c + (7 a - 7 C ) sin 2 9 n , 
j'=± 



(9) 



(10) 



for the total transition probabilities from the states |±,n). 

Although the transitions occur between states of different 
degrees of entanglement, the total transition probabilities are 
determined only by the amplitudes of the states of the upper 



doublet. The states of the lower doublet do not become in- 
volved. The properties of the total probabilities are essen- 
tially similar to those discussed above for the transitions from 
the n = 1 doublet to the ground state. Indeed, for n = 1 
the probabilities ( TTOb reduce to Eq. Q, and unequal damp- 
ing rates, 7 a 7^ j c , ensures the dependence of the transition 
probabilities on the degree of entanglement. 

Thus, transition rates of a combined quantum system de- 
pend strongly on how the subsystems decay rather than on 
the nature of states involved. The presence of the maximally 
entangled states in a JC cell could be observed in principle 
by measuring the transition rates between energy levels of the 
system, subject to 7 a 7^ j c . 



C. Absorption spectrum of a weak probe beam 

There remains the question how these properties of the tran- 
sition rates might be exhibited experimentally. Probe-beam 
spectroscopy methods should be able to test these properties. 
When a JC system is irradiated by a weak probe beam with a 
frequency that is close to resonance, the rate the probe beam 
is absorbed is proportional to the transition rates between the 
energy states of the system. 

We consider net absorption by the system of radiation from 
a tunable beam probing the system. The probe beam intensity 
is assumed to be sufficiently weak such that the field is treated 
to only first order in its amplitude so that it does not apprecia- 
bly perturb the system. The absorption spectrum is given by 
the ima gina ry 
system Il8l420ll 

7 r 



art of the linear susceptibility Xo (w p ) of the 



Im 



X^K) 



(oj in -w p ) + 7 2 



(11) 



where oj p is the frequency of a probe beam, T^n is the total 
transition rate from the state i (i = ±) of the manifold n, 7 de- 
scribes width of a given transition. In writing expression (fTTT i. 
we have used the fact that transitions from the excitation states 
to states of the manifold below might not be purely radia- 
tive HUT]], i.e. 7 7^ T^n/2, and the system was in the n — 1 
manifold before the probe was applied. To ensure that the 
spectral lines are well resolved, we assume that the transitions 
do not overlap. This is achieved assuming that the coupling 
strength g is much larger than the transition rates r^. 

Figure [2] shows absorption spectra of a probe beam tuned 
in vicinity of the transition frequencies from the n = 1 dou- 
blet to the ground state of a single JC cell. We see that as 
long as A 7^ 0, the spectrum of the single cell is markedly 
asymmetric. This feature is associated with the fact that at 
A 7^ the transitions from the n = 1 doublet to the ground 
state occur with different rates, r + i 7^ T_.i. The spectrum 
becomes symmetric at resonance, where A = 0. The symmet- 
ric spectrum is associated with the fact that at the resonance, 
r + .i = r_,i. As predicted by Eq. (0, the equality of the 
total transition rates takes place only for the maximally entan- 
gled states. Therefore, the symmetry of the spectrum can be 
regarded as an indication of the presence of maximally entan- 
gled states. 
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FIG. 2: (Color online) The absorption spectra Im £ 
probe beam monitoring a single JC cell for 7a/<? = 



of a 

1 = 0.05, 7c /<7 = 
0.02, j/g = 0.01, and different detunings: A/g = (red line (iii)), 
A/g = 2 (black line (ii)), and A/g = 10 (blue line (i)). 



III. TWO MUTUALLY COUPLED JC CELLS 



is degenerate 



(Dl 



,JC(0) 



,JC(1) 



(14) 



and the spectrum corresponds to what is expected for two in- 
dependent JC systems, with e being irrelevant. 

On the other hand, for g = = A, the v = spectrum 
is characterized by uj^ = 0, whereas the v = 1 is described 

by uril^ | s= o=a, which is the doubly-degenerate lo c or ui c + en 
values as expected for coupled harmonic oscillators. There- 
fore, a large inter-cavity hopping rate will push the coupled 
JC cells to behaving closely like coupled harmonic oscillators 
with atom-cavity perturbations. For small k 



(15) 



with shift ±g due to vacuum Rabi splitting, ek/2 the normal- 
mode splitting due to inter-cavity coupling and en 2 /8g a fre- 
quency pulling effect similar to the ac Stark shift. 



We now consider two neighboring JC cells coupled via 
overlapping evanescent waves of the cavity modes. The cou- 
pling results in an coherent hopping rate k between cavities. 
We designate Sj and rjj as the field annihilation operator and 
atomic electron energy lowering operator for the i' 
and the double JC cell Hamiltonian is HH-O 



th JC cell, 



H = H 



JC 



k ( a\a,2 + &1&2 



(12) 



for v denoting the total number of quanta shared in the two- 
cavity system. 

The space of the system for any v is spanned by prod- 
uct states {|ni, ci, ri2, C2)} with Tij a label for the number 
of photons in the i th cavity and c; the label for whether the 
state is in the excited (|e)) or ground (\g)) state. For ex- 
ample, the ground state corresponds to v = with just 
one dimension hence contains one basis element { | O.gOc?) }. 
The single-excitation (v = 1) basis comprises four states 
{\0e0g) ,\lg0g) ,\0glg) ,\0g0e)}. 

In the following, we limit ourselves to the subspace v = 1 
that only a single excitation is present in the system. The v = 
1 case is especially interesting as it corresponds to four qubits: 
two two-level atoms and two field modes with superpositions 
of vacuum and single-photon states. Also the v = 1 case can 
be solved in closed form yielding simple expressions Q. 



A. Entangled four-qubit states 

A diagonalization of the Hamiltonians for v = 0, 1 
yields the energy spectra = —A, for the ground state, 
and 



i(A + e K ) 2 (13) 



for the single-excitation, with e = ± and e = ±. In the case 
of independent cells, k = 0, and then the upper spectral value 
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FIG. 3: (Color online) The dependence of the energy difference 
8u}^ — Wei — on the detuning A (in arbitrary units) of 
the v = 1 states, | + , — ) (black solid line), |— , — ) (red long dashed 
line), |— , +) (green short dashed line), |+, +) (blue dotted line), for 
two coupled JC cells with (a) g = 0, n — 2, and (b) g — 1, k = 2. 

Figure[3]shows the variation of the energies lo^ of the sin- 
gle excitation v = 1 states with the detuning A around the 
unperturbed energy cj c . In the absence of the coupling g, there 
are two crossing points at A = ±k. The coupling lifts the 
degeneracy resulting in four entangled states. The states are 
eigenstates of the Hamiltonian ( fT2l ) corresponding to the en- 
= ergies u>il , and are of the form [0] 

|e, e) = u c , e (\lg0g)-e |0 5 l 5 ))+ We , £ (|0e0 5 )- e|0 5 0e}) , 

(16) 

where 



—TV + £ 



+ £v /TT72 



2 + 2 r e -ey/l + r* 



2 + 2 [r e -e^/T 



and 



r e = (A + os)/(2fl). 



(17) 



(18) 



The double-JC states ( TToT i are in the W state class correspond- 
ing to a superposition state of single excitations amongst each 
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of the two atoms and two modes. These states are maximally 
entangled only for w e e = w e . e = which requires that 
r t = 0. The condition that r t = is met only if A = en. 
Radiative properties of two coupled JC cells are studied in the 
next section. 



B. Transition rates and their collective properties 

We can compute the transition probabilities lfl3l - [l5ll 

T e , s = 7a |<e, s\ (at + <r+) |0) | 2 + 7c | (e, e\(a\ +4) |0) | 2 (19) 

from the single excitation states |e,e) to the ground 
state |0g0g) = |0), which yields 

Te.e = (1 - ef ( 7a \w t , £ \ 2 + 7c K £ | 2 ) . (20) 

The dependence of the transition rates on e signals the exis- 
tence of superradiant (e = — ) and metastable non-radiative 
(e = +) states in the system lfl6i [l7il . Note that the equal- 
ity of the superradiant rates, = r_._, is a signature of 
the maximal entanglement present in the corresponding states. 
However, this signature is not present for the other two states 
|+,e), as these do not imprint a signature on the radiation 
properties of the system. 

Notice that these collective radiative properties of the sys- 
tem are independent of the ratio 7 a /7 c and of the nature of the 
states. These decay processes are due to by losses of the cou- 
pled cavity modes as well as by atomic spontaneous emission. 
This latter feature is surprising because atoms are not directly 
coupled to each other, and they behave independent of the ra- 
tio n/g. The existence of the metastable states implies that the 
double JC system effectively behaves as a single collective JC 
cell composed of a doublet consisting of four-qubit entangled 
states |e, — ). 

The metastable states could be made radiatively active by 
breaking the symmetry between the atoms and/or the cav- 
ity modes of different JC cells, for example, by allowing the 
atoms and the cavity modes to decay at different rates. It is 
easy to show that when the atoms are damped with differ- 
ent rates, say 7 a i and 7^, and the cavity losses occur with 
rates 7 c i and 7 C 2, respectively, then the transitions from the 
states |e,e) occur with probabilities 

r e , e = (VT^l - ^\/l^l) 2 \We.e\ 2 + (V^tl - ^y^) 2 |w e , e | 2 • 

(21) 

Clearly, the transition probabilities are different from zero ir- 
respective of e. However, the states retain their collective char- 
acter with the states e = — still behaving as a superradiant and 
the e = + states now behaving as a subradiant state. 

There are some similarities in the properties of the transi- 
tion probabilities of the double- and single- JC cell systems, in 
particular, when the atomic spontaneous emission and cavity 
losses in a JC cell occur with the same rate. For example, in 
the case of 7 a i = 7c i and 7a 2 = 7c 2, the probabilities (ETb 
become independent of the amplitudes of the states, which is 
the same property encountered for J?±i of the single cell. 



However, there are interesting differences. In the double JC 
system the condition 7 a i 7^ 7 C ; is necessary but not sufficient 
for the dependence of the transition rates on the amplitudes of 
the states. There is also a rather subtle condition of the relation 
between the damping rates of different cells to be satisfied. 
Even though the spontaneous emission and cavity losses in a 
given JC cell occur with different rates, 7 a j 7^ 7^, the transi- 
tion probabilities still could be independent of the amplitudes 
of the states. This happens when 7 a i = 7 C 2 and 7 a 2 = j c i, for 
which Eq. (|2TT> reduces to 

1 2 

r e , e = - {y/^hi - eVTci) ■ (22) 

Evidently, the transition rates are independent of the am- 
plitudes of the states. The quantitative reason for this is that, 
even if the cells are distinguished by different damping rates, 
they are directly coupled to each other through the coupling k. 
This coupling creates entangled states between the cells that in 
the case of 7 a i = 7 C 2 makes indistinguishable through which 
channel, 7a i or 7 C 2, a photon leaving the system was emitted. 
This property is characteristic of multi-cell JC systems and 
does not exist in a single JC cell. Therefore, under the condi- 
tion that all the damping rates are different, equality between 
the superradiant (subradiant) rates signals maximal entangle- 
ment in the superradiant (subradiant) states, |— ,e) (|+,e)), 
which generalizes the results of the single JC cell. 

C. Absorption spectra 

We now consider absorption spectra of a weak radiation 
monitoring the transitions from the single-excitation (y = 1) 
states to the ground state |0) of two identical JC cells. To 
ensure that the probe field couples exclusively to the single 
excitation states, we assume that the transition do not overlap. 
This is achieved assuming that the coupling strengths g and k 
are much larger than the transition rates T m . 
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FIG. 4: (Color online) The absorption spectra Im ^Xt/' (<^ P )J of a 

probe beam monitoring two coupled identical JC cells, for 7a/g = 
0.05, 7 c /g = 0.02, j/g = 0.01, n/g = 2, and different detunings: 
A/g = 1 (red line (iii)), A/g = 2 (black line (ii)), and A/g = 10 
(blue line (i)). 

Figure|4]shows absorption spectra for 7 a 7^ j c and different 
detunings A. We see that as long as A 7^ ±«, the spectrum is 
always asymmetric. The spectrum becomes symmetric when 
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A = «. In this case the states |±, — ) are maximally entangled 
states. Thus, similar to the case of a single JC cell, the sym- 
metry of the spectrum can be regarded as an indication of the 
presence of maximally entangled states. Note that the sym- 
metric spectrum is observed at non-zero detunings, A = ±k, 
which is in contrast to the properties of the spectrum of a sin- 
gle cell, where the symmetric spectrum is observed only at 
resonance, where A = 0. 

(b) 
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FIG. 5: (Color online) The absorption spectra Im [xo^( w p)] °f a 

probe beam monitoring the system of two nonidentical IC cells for 

(7ai/<7,7a2/ff) = (0.01,0.2), (7ci/3,W.g) = (0.2,0.05), n/g = 
4, y/g = 0.05, and different A: (a) A/g = -4 and (b) A/g = 4. 

Figure [5] shows the absorption spectrum for a more general 
case of unequal damping rates of the atoms (y a \ ^ 7^) and of 
the cavity modes (7 c i 7^ 7^). In this case, the spectrum com- 
prises four peaks corresponding to the transition rates of the 
four eigenstates |e, e). Two of the peaks are high (correspond- 
ing to the superradiant states | — , e)), and two are short (cor- 
responding to the subradiant states |+, e)). According to our 
predictions, whenever the high (short) peaks have the same 
height and are located in opposite sides around the cavity fre- 
quency, this condition corresponds to a signature of the under- 
lying maximal entanglement of the superradiant (subradiant) 
states. The joint y a \ = y a 2 and y c \ = y C 2 case is special be- 
cause, in this case, the subradiant states are optically inactive 
and only the peaks corresponding to the superradiant states are 
present in the spectrum. Therefore, the symmetry of the spec- 
trum around the cavity frequency increases as the superradiant 
states \—,e) become increasingly entangled. Note that, in or- 
der to resolve the peaks, the linewidth of the probe beam has 
to be sufficiently narrow compared with the cavity and atomic 
damping rates. 



D. Atoms damped by independent reservoirs 

In the above analysis we have assumed that the atoms and 
the cavity modes are damped by common reservoirs. This as- 
sumption is justified if the cells are quite close to each other 
where the spatial variation of the field can be ignored. This as- 
sumption may not always be true in practice as the cells could 
be separated by a large distance. In this case the assumptions 
that the atoms are coupled to the same reservoir might not be 
true. For this reason, we consider now a situation illustrated 
in Fig. [6] in which atoms located in distant cells are damped 
by independent reservoirs. 

In the case of independent reservoirs for the atoms, the tran- 
sition probabilities from the single excitation states \e,e) to 




FIG. 6: (Color online) Two distant cells with the cavity modes over- 
lapping in the shadow area. This area can also be treated as a com- 
mon reservoir for the cavity modes to which the modes are damped 
with rate y a . The atoms located in the distance cavities are damped 
by independent reservoirs (green areas) with rate y a . 



the ground state |0) are given by 

r e , e = 7 a (|( £ ,e|a+|0)| 2 + |(6, £ | ( 7 2 + |0)| 2 ) 



+ 7c|( e , £ |(a t 1 +4)|0)| 2 , 



(23) 



which yields 

r,. 



27aKe| 2 + (l-e) 2 7cKe| 2 . ( 24 > 



We see that, in the case of independent reservoirs for the 
atoms, the transition probabilities T e e all become different 
from zero. There are no trapping (subradiant) states that 
would form a subspace of decoherence-free states. However, 
in the limit of |u> ei£ | 2 = 0, the states e = + would form a 
decoherence free subspace. 

A close look at Eq. ( TT7T > reveals that the amplitudes \w e> - 1 2 
vanish in the limit of r e 3> 1. The condition of r e 3> 1 is met 
when (A + e«;) 3> g which shows that either the atomic transi- 
tion frequencies differ significantly from the cavity frequency, 
A ^> 0, or the cells are strongly coupled to each other, n ^> 0. 
Under this condition, 



and 



K+I 2 = K-| 2 



0, 



(25) 



(26) 



As a consequence the states of the system in Eq. ( TToT i reduce to 
1 



| e ,+) = _(| e5 )- e | 5 e))|00) 



1 



V2 



(|10)-e|01))| 55 ) 



(27) 



with the corresponding transition probabilities 

r e ,+ =7a, T £ ,_ = i(l-e) 2 7c . (28) 

Clearly, the transition probability = 0, so that the 

state 1+, — ) is a decoherence free state. 



We stress that, in contrast to the states ([Tol l damped by 
a common reservoir, the states (l27t are product states with 
atom-cavity decoupling due to having independent reservoirs 
and strong inter-cell coupling. Effectively, the system behaves 
as a two-qubit system. Our prediction agrees with the results 
of de Ponte et al. lfl2ll . who considered a system of coupled 
resonators interacting with independent reservoirs and found 
that the resonators must be strongly coupled to each other for 
the accomplishment of the condition leading to the decoher- 
ence free states. 



ric eigenstates. These are 



1 



\<h) = jjr P+\0g0g---0e)-\0g0g---lg) 



— \0e0g---0g) + \lg0g---0g) 



1 

\0g0g---lg) + \0g---lg0g)], 



\<h) =Jf[p+{ IO5O5 <></••• OeOg) 



IV. GENERALIZATION TO AN ARBITRARY NUMBER OF 
COUPLED JC CELLS 



Although the focus of this work is on two coupled JC cells, 
we can solve the problem more generally. In this section we 
show how to solve elegantly the case of N mutually con- 
nected JC cells with the Hamiltonian 



H 



N 



N 
i=l 



He 



(29) 



The Hamiltonian describes a system composed of N identi- 
cally coupled cells. While this description does not include 
such features as the direct coupling between the atoms and 
the spatial variation of the field modes, it is nevertheless of 
interest, as the simplest model of a group of collective behav- 
ing cells. This case is somewhat artificial for large N but does 
yield an instructive generalization of the formalism in addition 
to providing some connection with condensed-matter studies 
of coupled JC systems @]. 



A. Energy spectrum 

Restricting to the single-excitation, v = 1, the space of the 
system is spanned by 2N product states 



{|0 5 0s--.0e),|0fl0ff.--l 5 ) > . 
\0eOg---0g),\lg0g---0g)}. 



(30) 



In the basis of these states the Hamiltonian ( 1291 can be repre- 
sented by a 2N x 27V matrix. Despite a large size, the matrix 
can be directly diagonalized and results in TV-cell superposi- 
tion states (eigenstates). We find that among the 2N states, 
one can distinguish two sets of N — 1 degenerate antisymmet- 



•>N-l 



1 



p + (|0 5 5 ---0e)-|0 5 0e...0<?) 



-\0g0g---lg) + \0glg---0g) 
with corresponding eigenenergies 



(31) 



uji — ■ ■ ■ — lon-i 



= Wc-^(A-k) + Ji(A- K ) 2 + <? 2 , (32) 



where 



M = M 



p± = —r ± \/ r 2 + 1, 



2, Af 2 = J2pl + 2, 



A — K 



ig 



(33) 



The other set of N — 1 degenerate antisymmetric eigen- 
states is 

I'M - j^r [p. \0g0g ■ • • Oe) - |0 5 5 • • • lg) 



\<f>N + i) = jr 



— lOeOff • - - Off) + jlffOff • - - Op) 
P+ 

r p_(|0«70 5 ---0 e )-|0 5 ---0e0 5 )) 



V2N-2 



1 



p_(|0 5 0ff---0e)-|0 5 0e---0 5 )) 



-\GgOg---lg) + \Qglg---Qg) 
with eigenenergies 



(34) 



Wjv = • • • — UJ2N-2 



= w c -i(A-«)- Ji(A-ft) 2 + .9 2 , (35) 



where 



A/3 = Jp 2 _ + 1- + 2, A4 = 1 
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The remaining two eigenstates of the system are fully sym- 
metric superposition states 



1 



1 



l&t) = tt — l°30ff ' ' ' Oe) + \0g0g ■ ■ ■ OeOg) +■■■ 
+ \0e0g ■■■0g}) + \0g0g ■■■lg) + \0g0g ■ ■ ■ IgOg) 



+ --- + \lg0g---0g)\, 
with eigenenergies 

w± = w c - T [A + (J\T-1)k] 



(37) 



A /i[A + (7V-l) K ] 2 +. 9 2 , 



respectively, where 




and 



± 



A + (/V- 1)k 
2^ ' 



(38) 



(39) 



(40) 



We note from Eqs. <|32), ([35} and (Q8J that for g = 0, en- 
ergies of the antisymmetric states cross at A = n, whereas 
energies of the two fully symmetric states cross at A = 
— (N — 1)k. When g ^ 0, avoiding crossings occur at A = k 
for the antisymmetric states, and at A = — (N — 1)k for the 
symmetric states. If and only if A = — (N — l)«is true, do the 
eigenstates \4>±) become W-type maximally entangled states. 




-15 -10 -5 5 10 15 
A [arb. units] 



-15 -10 -5 5 10 15 
A [arb. units] 



7: (Color online) The dependence of the energy difference 
= uii — w c on the detuning A (in arbitrary units) of the v = 1 



FIG 

eigenstates of TV = 10 coupled IC cells for k = 1 and (a) g = 0, 
(b) g = 1. The green short dashed line is for the set of antisymmet- 
ric states \(f>i) (i = 1,2, . . . , N — 1), the red long dashed line is for 
the set of antisymmetric states \(j>j) (J = N, N + 1, . . . , 2N - 2), 
and the blue dotted and black solid lines are for the symmetric states 
\<f)+) and \4>-), respectively. 

Figure|7]shows the dependence of the energies of the eigen- 
states \cj>±) on the detuning A. We see two crossing points 
when .9 = and the avoided crossing between the states when 
g 7^ 0. It clearly illustrates the appearance of two separate 
groups of degenerate antisymmetric states with the minimum 
energy separation 2g at A = k, and two symmetric states with 



the minimum energy separation 2g at A = — (N — 1)k. Per- 
haps the most interesting aspects of the N cell system is that 
the energies of the antisymmetric states are independent of N. 
Adding more cells increases the number of the antisymmetric 
states but does not affect their energies. Here the antisymme- 
try is with respect to a permutation of the two entangled JC 
cells, not of any pair of 'qubits', and it is exact only in the 
limit of maximal entanglement \p±\ = 1. For the specific 
case N = 2, the antisymmetric states |(/>i,2) correspond to 
the (subradiant) states |+, i) introduced in Eq. (TToT l. whereas 
the symmetric states \ <fi±) correspond to the superradiant ones 
| — , i). Note also that N = 2 is the only case in which the an- 
tisymmetric and symmetric states have entanglement between 
the same number of 'qubits', four in this case. 



B. Transition rates 

Having derived the explicit forms of the eigenstates \<j>i) of 
the single-excitation sector v = 1, we now turn to calculate 
transition rates between the eigenstates \<pj) and the ground 
state |0) = \Qg ■ ■ ■ 0g). We consider separately two cases. In 
the first, we assume that the atoms are coupled to a common 
reservoir. In the other case, we assume that the atoms are 
coupled to independent reservoirs. 

When the atoms are coupled to a common reservoir, the 
transitions occur with rates 



r, = 



N 

i=i 



(41) 



Applying Eqs. PTT i. (|34| | and (l37l i. after straightforward cal- 
culations we obtain 



i 

r, 




— \ 2 

dN \ it i — /— v 

—J +(Vt^-V7 c1 ) 



7ai) 2 + (\/7cI - \/7 i)' 



+ (VTfew- V7 c i)' 



j = N + 1,...,2N -2, 



± 



1 



For the case of identical cells, 

7al = • • ■ = laN = 7a 



and 



7ci 



7cJV = 7c 



(42) 



(43) 



(44) 
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it follows that the transition rates from all of the antisymmetric 
states are zero 



N 



2N-2 



0. 



(45) 



Only the fully symmetric states \4>±) 
are damped at rates 



are optically active and 



r+ = 



TV 



(I+P'l 



75T (7a+P±7c), 



(46) 



in which we see the characteristic N time enhancement, the 
superradiant behaviour [fl6l l22ll . 

From Eqs. (|43T > and (|4o*T i it is apparent that, just as in the 
N = 2 case, the 2N — 2 antisymmetric states are all op- 
tically inactive, whereas the two symmetric ones are super- 
radiant. The distinction is quite convenient because higher- 
dimensional entanglement is present in the latter states, so a 
symmetric absorption spectrum would reveal maximal 2N- 
partite W-like entanglement shared among all the JC cells. 

As the TV-cell system radiates only from the fully symmet- 
ric states, the Hilbert space of the system is confined to the to- 
tally symmetric subspace of two states \4>±). Given that there 
are only two non-zero transition rates from the excited states 
to the ground state, the absorption spectrum of a weak probe 
beam monitoring the system is expected to be composed of 
two peaks. This feature is seen in Fig. [8] which displays the 
absorption spectrum for coupled N = 10 cells. 



20 

^ 15 
1 10 
5 




(i) 



(iii) 



(ii) 



•10 



-5 

(w c -w p )/g 
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FIG. 8: (Color online) The absorption spectra Im ^Xo^C^p)] of 
coupled iV = 10 cells plotted as a function of (cj c — u v )/g for 
7a/ff — 0.05, 7c/.g = 0.02, n/g = 2, and different detunings 
A: A/g — 1 (green line (iii)), A/g = 10 (black line (ii)), and 
A/g = -18 (red line (i)). 

As expected, there are two peaks corresponding to transi- 
tions from the states \<fi±) to the ground state |0). The rela- 
tive peak amplitudes are a function of A, and if and only if 
A = — (N —1)k, the amplitudes are equal and symmetrically 
located about 



uj c -u> p = ±g. 



(47) 



When A = — (N — 1)k, the states \<f>±) reduce to maximally 
entangled W states. Hence, a symmetric absorption spectrum 
is a signature of a maximally-entangled W states with 2N de- 
grees of freedom. 



We close this section by evaluating the transition rates for 
the case when the atoms are damped by independent reser- 
voirs. The total transition rate due to the interaction of 
the atoms with independent reservoirs and the cavity modes 
damped by a common reservoir is defined by 



N 

E^- 



E 



TtJ(&|aJ|o> 



(48) 



When we make use of Eqs. (1311) . ( 134b and d37l i in Eq. (|48}, we 
readily obtain 



r, 

i 

r« = 
3 



1 
1 

TV? 




A/? 



Ml 

N 4 



p\ (7al +7a0 + (V7 

,N-1, 

(p_7al + -!pr) + {Vi^-Vici) 

P- (7ai + 7<w) + (V>7 - \/7 c i) 
l,...,2iV-2, 



1 

M 



1 ^ (» 

F± i=l \i=l 



(49) 



It is easy to see that the transition rates (149) differ significantly 
from that we encountered in Eqs. (l42l i for the decay of the 
atoms in a common reservoir. Notice that all of the transition 
rates are now different from zero, and none of the rates can 
be reduced to zero, indicating no subradiance. This principle 
also applies to the special case of identical cells, 



and 



7al = ' ' ' = 7aW = 7a 



7cl = ' ' ' = 7cW = 7c 



(50) 



(51) 



However, despite the fact that the transition rates (149) are 
different from zero, we find that in the case of identical cells 
and under a strong coupling between the cells (k ^> g), the 
transition rates reduce to 



Ti 



7a: 

0, 

7a, 



2, 



N, 



,N-1, 
..,2N-2, 



1 



P'l) 



(7a+P±iV7c 



(52) 



Clearly, in the limit of a strong coupling the states |<^) (i = 
2, . . . , N — 1) do not decay, so they form a decoherence-free 
subspace. The remaining states decay with nonzero rates. The 
states \ (j>i) , \4>n) and \<pj) decay with the rate equal to the sin- 
gle atom decay rate y d , whereas the rates T± show the charac- 
teristic N times enhancement (superradiance), but only with 
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respect to the cavity damping. The contribution to r± from 
the atomic part, proportional to j. d , is independent of N, This 
is what one could expect, since the atoms are damped by in- 
dependent reservoirs. This is also consistent with the results 
of de Ponte etal. lfl2ll . 

V. CONCLUSIONS 

We have studied a system of coupled Jaynes- 
Cummings (JC) cells with coupling created by overlapping 
evanescent cavity fields or by tunneling of photons between 
the cells. This overlap or tunneling determines the coherent 
hopping rate between the pair of cavities. This coherent 
hopping term along with the strength of atom-field coupling 
within each cavity determines the unitary dynamics of the 
closed coupled system. Our analysis also includes incoherent 
cavity damping and atomic spontaneous emission rates. 
These incoherent processes are unavoidable in experiments 
but are sometimes not considered in studies of condensed- 
matter types of properties of such systems. We included 
these incoherent terms for completeness and found that 
ensuring differences in incoherent rates can be quite valuable 
in detecting entanglement of these coupled atom-cavity 
systems. 

The principal advantages of our treatment are that it allows 
to study the relation between the radiative properties of the 
system and entanglement. For example, in the case of the sin- 
gle JC cell we have observed that the ratio of the transition 
rates for the two lowest-doublet states directly reveal the de- 
gree of entanglement between the atom and cavity mode pro- 
vided that the spontaneous emission rate 7 a is not equal to the 
cavity loss rate 7 C . The degree of entanglement depends on the 
values of atom-cavity coupling strength g in the JC Hamilto- 
nian. 

We have also studied transition rates from the second dou- 
blet to the first doublet for the single JC model. In this case 
we have transitions from entangled n = 2 state to entangled 



n = 1 states, in contrast to the n = 1 to n = transi- 
tion from an entangled first-doublet state to a ground state, 
which is an atom-cavity product state. We have found that 
the rates depend explicitly on properties of states in both dou- 
blets, namely that it depends explicitly on both 9 n and 6 n -\, 
but adding the transitions rates together as in Eq. (TTOt only de- 
pends on 6 n . Thus the transition rates directly reveal the de- 
gree of entanglement of the n th doublet without being mixed 
up with the degree of entanglement in the lower doublet. Di- 
rect measurement of entanglement of the states is possible. 

These results on entanglement in the single JC cell are inter- 
esting and provide a foundation to study coupled JC cells. For 
coupled JC cells we focus only on the v = 1 case. We have 
derived the explicit analytical forms of the transition rates be- 
tween the v = 1 quadruplet of states for the double JC cell to 
the (product) ground state. There are four rates in this case, 
and each of these rates conveys information about the degree 
of entanglement in the state including whether the states are 
maximally entangled. When loss rates from cells (both spon- 
taneous emission and cavity loss) are equal, information about 
the degree of entanglement is not obtained by our method but 
is obtained for unequal rates. 

Finally, we have shown how to extend our two JC cell case 
to an arbitrary number N of the JC cells such that each cell is 
coupled to all other cells. We have established conditions for 
maximal entanglement and calculate the susceptibility. This 
N cell case is likely only of theoretical interest as it is hard 
to imagine how such a system would be configured in a real 
experiment, but this general extension to N cells is interesting 
nonetheless. 
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